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SUMMARY 


A new class of nonideal solutions is defined by constructing a function 
to represent the composition dependence of thermodynamic properties for mem- 
bers of the class, and some properties of these solutions are studied. The 
constructed function has several useful features: (1) Its parameters occur 

linearly. (2) It contains a logarithmic singularity in the dilute solution 
region and contains ideal solutions and regular solutions as special cases. 
(3) It is applicable to N-ary systems and reduces to M-ary systems (M <_ N) 
in a form-invariant manner. The ability of the function to describe data 
for a real system is tested in a companion report (NASA TP-1930). 

INTRODUCTION 

The effective application of classical thermodynamics to complex physi- 
cal problems demands two things. First, it requires an understanding of the 
foundations of thermodynamics in order to appreciate both its capabilities 
and its inherent limitations. Second, it demands convenient and accurate 
mathematical representations of thermodynamic data in order to facilitate 
their use in computations. On the one hand, theoretical studies of thermody- 
namics are gradually and successfully exposing its foundations (ref. 1). On 
the other hand, theoretical thermodynamic studies have had little real suc- 
cess in guiding us to an accurate representation of thermodynamic data. The 
many different functions which have been, and are being, used to represent 
the composition dependence of thermodynamic properties are ample evidence of 
this. Each of these many functions, some wholly empirical and others with a 
theoretical base, has been used to give an adequate representation of data 
only for a limited number of systems or for a restricted composition range. 
None of them seems to have the flexibility to deal with a large variety of 
systems over their entire composition range. It would be a great computa- 
tional convenience to have available a relatively simple "universal" function 
of composition variables which is capable of adequately representing the data 
for all systems over the complete range of compositions and which can be used 
for interpolation and extrapolation of the data. As nice as this would be, 
we have no assurances either that such a function exists or, if it exists, 
that we can find it. In fact, given the diversity of physical systems, it 
is more likely that such a universal function cannot be found. However, it 
does seem reasonable to expect that we can construct a function which is, in 
this respect, an improvement over those already being used. 

In this report I shall first briefly survey some of the composition 
functions which have been used to represent thermodynamic data. Then I will 
construct a new function which I propose for the representation of data. 
Finally, I shall examine some of the properties of the class of solutions 
which can be represented by this function. The application of the new 
representation to the experimental data for a few different systems is dis- 
cussed in a companion report (ref. 2). 

In making a survey of the functions which have been used for the compo- 
sition dependence of thermodynamic properties one quickly realizes two 
things. First, most of the functions tend to be expressed in terms of poly- 
nomials; second, by far the greatest variety of functions has been applied 
to binary systems, with ternary systems being in a distinct second place. 

This is a natural consequence of the need to deal with the idiosyncrasies of 
actual data for a particular system and the fact that experimental data have 



been pretty much confined to binary systems with a small sprinkling of tern- 
ary systems. Marsh (ref. 3) has observed that many of the approximations to 
the thermodynamic excess functions of binary systems can be regarded as 
special cases of a rational function of two polynomials. Thus if aG is 
the excess Gibbs free energy per mole, T the absolute temperature, R the 
universal gas constant, and xi,x2 mole fractions, then he considered 
functions of the form 


aG 

Xj^RT 


m k 

Z) M X 1 " 

k=0 K 1 c 

n 

Z) B.-fx, - x ? ) J 
j=0 J 1 * 


[n,m], B q = 1 


( 1 ) 


Of course, for binary systems X 1 + *2 = 1 and equation (1) could be written 
as a function of x^ or X 2 alone. Approximations of the type (1) can 
clearly be regarded as Pad§ approximants (ref. 4) of numerator degree m and 
denominator degree n for aG/xix 2 RT, and this gives us a convenient classi- 
fication scheme for discussing tne many expressions used for binary systems. 
The approximants [0,m] and [n,0] have been discussed by Van Ness (ref. 5). 
Otterstedt and Missen (ref. 6) refer to the approximant [0,m] as the general 
Margules representation and refer to approximants of the type [n,0] as the 
general van Larr representation because [1,0] corresponds to the original van 
Laar form (ref. 7). Myers and Scott (ref. 8) used the approximants [l,m], 
m <_ 4, which might be regarded as modified Margules representations although 
they could equally well be viewed as generalizations of the van Laar form. 

Van Ness and coworkers suggested three equations which they regarded as mod- 
ifications of the Margules form. Abbott and Van Ness (ref. 9) suggested the 
forms [1,2] and [2,3], while Shatas, Abbott, and Van Ness (ref. 10) proposed 
the form [3,4] - all are of the type [m - 1, m]. The first of these is a 
special case of the functions used by Myers and Scott. Marsh expressed the 
second of these as a ratio of two quadratics but, from the Abbott and Van 
Ness version of the formula, it is clear that the numerator must be a cubic 
polynomial and therefore Marsh's formulas relating the coefficients in equa- 
tion (1) to those used by Abbott and Van Ness may be incorrect. Not only 
can approximants of the type [n,0] and [l,m] be regarded as generalizations 
of the van Laar form [1,0], but also there exists a third type which can be 
viewed logically in the same way. These are the approximants [n, n - 1] 
since they too contain [1,0] as a special case. These approximants also in- 
clude as special cases the "q-equations" of Wohl (ref. 11) for binary sys- 
tems. His (n + l)-suffix q-equation corresponds to the approximant [n, n - 1] 
if the coefficients of the denominator polynomial are of the form 



where 



is a binomial coefficient and the 


on the right is a parameter expressible in terms of Wohl's parameters q^ 


B 


and q£. Thus B = ( q^ - q^ ) / ( q^ + ^ ’ Wo ^ actually considered only the 

cases n = 1, 2, 3, 4. There is still another approximation of binary 
properties which should be mentioned. It is the square root function of 
Otterstedt and Missen (ref. 6) and is not of Pad£ approximant type. 
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Otterstedt and Missen use = 1 and Bandreth (ref. 12) uses n^ = 2. 

The functions which have been mentioned to this point are only appli- 
cable to binary systems. If we wish to consider systems with N constitu- 
ents (N-ary systems), then we must use other forms which, of course, will 
apply to the binary forms for N = 2. The most obvious expression one could 
write would be an analog to a truncated Taylor expansion, and it could be 
thought of as a virial expansion. 
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Benedict, Johnson, Solomon, and Rubin (ref. 13) used an equation of this 
type; however, they restricted it to use only the quadratic, cubic, or 
quartic terms. Not all of the coefficients in equation (3) are independent 
since the coefficients must be symmetric in the indices. Also since aG is 
an excess function, then for all i A + A-j + A-j-j + A-p-j + ... = 0. Another 
form which can also be interpreted as a truncated Taylor series but with 
different term groupings is 
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This form of representation was used by Pelton and Bale (ref. 14) for a 
ternary system (N = 3). The coefficients in equation (4) must satisfy 
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for excess functions. A variant of equation (3) which has been given consid- 
erable attention is the approximation 
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(5) 


For the choice A-j = 0 this becomes the N-ary version of Wohl's (ref. 11) 
q-equations for ternary systems, which he wrote through quartic terms. 
Another well-known special case can be obtained from equation (5) if we 
choose Aj, A-jjk, and all higher-order coefficients as zero. Supple- 
menting these conditions with = 1 and A^j = w(l - 6 i j)/2, where 
fi nj = 4j& is the Kronecker delta, gives the excess Gibbs free energy for 
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the regular solutions of Hildebrand (ref. 15). If instead we pick A^j 
= w(i - j/2, then we obtain an N-ary generalization of the binary 

van Laar solutions, that is, those characterized by the approximant [1,0]. 
Perhaps the most extensive use of equation (5) has been by Pitzer and 
coworkers (ref. 16), who applied it to aqueous electrolyte solutions. Their 
excess function can be obtained from equation (5) by selecting A-j = 0 
for all i, B-j = 0 for i ^ 1, and equal to the molecular weight of 
water. The form of the first term in equation (5) also suggests an N-ary 
analog of the Pad£ approximants for binary systems, that is, the ratio of 
two functions of the kind shown in equation (3), although I know of no 
instance where such an approximation has been attempted other than the two 
special cases just mentioned. 

Each of the functions mentioned to this point has been firmly rooted in 
the polynomials. This is not true for the next example. 

r N / N \ N/N /N \ 

RT - A £ x i (jg A i j x j J + B £ x i B i J-'j/S C ik x kj (6) 

This function was given by Renon and Prausnitz (ref. 17), who actually 
wrote it with C-j^ = A-j^. It contains a number of functions, proposed 
by others, as special cases. If B = 0, it becomes the Wilson equation 
(ref. 18). It becomes the equation of He i 1 and Prausnitz (ref. 19) for 
B = -A = 1 and C -j ^ = Aj [< and the NRTL (nonrandom, two liquid) equation 
of Renon and Prausnitz (ref. 17) for A = 0 and B = 1. 

A final illustration of the kinds of functions which have been used to 
represent the composition dependence of thermodynamic data is one which 
might be called the "bootstrap" representation. It is not as specific a 
functional form as those I have already mentioned, and it could utilize any 
one of them. Although this representation could be written for an N-ary 
system, the notation becomes complex and so I shall write it only for a 
ternary system. 

aG ( X-^ , X 2 , X 3 ) = aG.^( x l ,x 2 ^ + A ^i 3 ( x i » X 3 ) + ^23^ x 2’ x 3^ + ^i23(*i’^2’ x 3^ (^) 

In this expression aG-,- j(x-j,x j) are any suitable binary functions and 
aGi23( x i,X9, x 3) is a ternary function. I call this a "bootstrap" repre- 
sentation Because it is possible, but not necessary, to determine the coef- 
ficients in the binary functions from binary data alone and subsequently the 
coefficients in the ternary function only from those ternary data for which 
x-j ^ 0 for all i. Shatas, Abbott, and Van Ness (ref. 10) used this kind of 
representation. For two of the binaries they chose a six-constant Margules 
representation from [0,5], and for the third binary they selected the [3,4] 
Pad£ approximant. The ternary function had the form of the cubic term in 
equation (3) multiplied by (x-^x^). For a quaternary system equation (7) 
would be replaced by the sum of a quaternary function and all possible binary 
and ternary combinations. 

It would probably be helpful to point out, in general terms, some of the 
difficulties connected with the use of the functions I have been discussing 
before attempting to construct an alternative. I have already mentioned that 
many of the functions lack the flexibility to handle the wide diversity of 
behavior encountered in real systems. One is then continually faced with the 
prospect of constructing new functions tailored to the system of interest or 
else using available functions in a piecewise manner as was done for alcohol- 
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chloromethane binaries by Moelwyn-Hughes and Missen (ref. 20) and Paraskevo- 
poulos and Missen (ref. 21). For most of these binary systems they found it 
necessary to use the approximant [0,2] for ln^G/x^RT) in the dilute alco- 
hol region and a three-constant Margules form elsewhere. For the binary 
methanol - carbon tetrachloride they used [0,2] in the dilute alcohol region 
and [0,1] elsewhere as approximations for In^G/x^RT) . Such piecewise 

fitting is a nuisance because one must ensure the continuity of the function 
and its first derivatives at the boundaries of the segments. The lack of a 
"universal" function is particularly disconcerting if one is attempting a 
Barker analysis (ref. 22) for the excess Gibbs free energy from vapor pres- 
sure data because one must make a commitment to a function in advance. The 
results are, to some extent, dependent on the function and its flexibility. 

A second source of difficulty is the fact that some of the functions I have 
mentioned are nonlinear in the parameters. This gives a person two options 
and neither is completely satisfactory. One can preselect some of the param- 
eters as did Myers and Scott (ref. 8) and Pitzer (ref. 16) and so reduce the 
equations to linear form by limiting their flexibility. Alternatively, one 
can face the unpleasantness of nonlinear parameter estimation and the possi- 
bility of nonunique solutions. A third problem arises because accurate ex- 
perimental data become increasingly difficult to obtain as one approaches 
the highly dilute solutions. But since the approximating equations are gen- 
erally unconstrained, they supply no guidance on the shape of the function 
in this region. It would be better if the function could be so constructed 
that it supplied some help in dealing with the dilute solution. Pitzer 
(ref. 16), for example, has done this in his work with aqueous electrolyte 
solutions by relying on Debye-Huckel theory as a guide. The final problem 
which I wish to mention is the possibility of introducing undesirable 
singularities into aG or its derivatives. This is built into the 
Otterstedt-Missen (ref. 6) formula and can arise in the functions with non- 
linear parameters by an inappropriate choice of parameters. The Otterstedt- 

- 1/2 

Missen function (2) with n^ = 1 produces an x^ singularity in the 

excess chemical potential for species 1 which dominates the In x-^ singular- 
ity from the ideal portion of the Gibbs free energy. 

In my opinion a suitable function for the representation of the compo- 
sition dependence of thermodynamic properties should exhibit six desirable 
features: (1) It should have a form directly applicable to N-ary systems. 

(2) It should reduce simply and obviously to lower order subsystems in a 
form-invariant manner. That is, the function for an M-ary system, M £ N, 
should have a form identical to that for an N-ary system. (3) It should 
contain ideal behavior as a special case because ideal behavior is often a 
good first approximation to real behavior. (4) Its behavior in very dilute 
solutions should be constrained to reflect the known behavior of these solu- 
tions. (5) Its parameters should appear linearly. Finally, (6) it should 
have sufficient flexibility to represent even difficult systems reasonably 
well . 

DEFINITION OF A CLASS OF NONIDEAL SOLUTIONS 

My desire to include ideal behavior as an integral part of the function 
rather than an add-on means that I must work with a thermodynamic potential 
rather than an excess function. Generally speaking, the most convenient 
experimental thermodynamic variables, other than composition variables, are 
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the intensive ones. These are the temperature and pressure in the case of a 
fluid. Consequently it is natural to work with the thermodynamic potential 
for which the moles of each constituent are the only extensive variables: 
the Gibbs free energy. In the following calculations only the composition 
variables will be written explicitly; the intensive variables will be sup- 
pressed. It should be kept in mind that all parameters are to be regarded 
as functions of the intensive variables. I shall begin by introducing nota- 
tion and by obtaining some results which will be useful subsequently. 

Consider a system of N constituents and let n-j, i = 1, 2, ..., N 
represent the moles of species i per unit mass. The Gibbs free energy per 
unit mass g is then a homogeneous function of degree 1 in the n-j. 

g(xn i ) = xg(n i ) (8) 

By definition the chemical potential for the i th species p-j is 

»f - 5 a 9 /sn i (9) 

and by Euler's theorem for homogeneous functions it satisfies 
N 

g = £ Vi (io) 

i=l 


It is a trivial consequence of equation (8) and the chain rule for differen- 
tiation that the chemical potential is a homogeneous function of degree zero 
in n-j. For x independent of n-j 


xy k ("i) = x [ag(n i )/3n |< J = jag(xn.)/an k J 

= £ [ag(xn i )/a(xnj)] [a( xn^/anj 


- x £ [ 3 9< xn ,)/3(*nj)]s jk 

- »rsg(>n j )/a(xn k )J . XM k (An j ) 
Hence we conclude that 

"k (n i> = " k Un,) 


( 11 ) 


and then from Euler's theorem and 


Van. = 


3 g/ an ^ 


an, 


N N 

£ n ap./an = 2J n. au./an, = 0 
j_l J K J J J K 


8 “j /8n k 


we obtain 


( 12 ) 
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The second equality is known as the Gibbs-Duhem equation. We now define the 
total moles per unit mass n and the mole fractions x-j. 


N 

n = E n • 
J=1 J 


x. = ni /n j 


(13) 


The choice x = n - ^- in equations (8) and (11) immediately shows that the 
chemical potential and the Gibbs free energy per mole G can be regarded as 
functions of the mole fractions. 


G = n 1 g(n i ) = g(n./n) = g(x i ) = G(x i ) = E 

j ^ 

p k = u k (n.) = v k (n./n) = y k ( Xi ) 


j 


(14) 


The last equality in the first member of equation (14) comes directly from 
equation (10). In principle, one can conclude from equation (14) that G(x^) 
and y|<(x-j) are homogeneous of degree 1 and 0, respectively, in the mole 
fractions because g(n-j) and y k (n.j) are suc ^ functions with respect to the 
moles. In practice these restrictions are both unnecessary and unwise be- 
cause they would force us to write unnecessarily complex functions. This is 

N 


a direct consequence of the fact that the mole fractions satisfy 
consider as an example the simple function 



n . 







Although the first two members are clearly homogeneous of degree 0, the last 
member is just as clearly homogeneous of degree 1. Over the domain of 


n 

definition E x. = 1, x. > 0, the last two members have the same value, but 
k-1 K J ~ 

N 

the values differ if extended outside this domain. The identity E x, = 1 


k=l 


can be used to interconvert mole fraction dependent functions which are homo- 
geneous functions of any degree and nonhomogeneous functions within the do- 
main of interest to thermodynamics, and so I shall disregard the homogeneous 
conditions. Because of equation (14) we can now express derivatives with re- 
spect to ni in terms of derivatives with respect to x-j. From equation (13) 
we find that 


ax^an^ = (6 ij . - x^/n (15) 

If we apply this and the chain rule for differentiation to equation (14), we 
obtain 
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}- ( 16 ) 


y i = 


1£_ 

9n i 


8[nG(Xj)] 



9x i 


G( Xj ) - (1 + *.)G( Xj ) 



n ' 1 Vk (x 1 ) 



where € ^ is a differentiation operator. The second member of equation 
(16), when combined with the second equality in equation (12), gives the 
Gibbs-Duhem equation written in mole fractions. 
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(17) 


We can now begin the actual construction of a function to represent the 
Gibbs free energy. Suppose we consider a system of N nonreacting species 
labeled by Latin indices i, j, k, ... = 1, 2, ..., N which can form two or 
more coe X isting phases labeled by Greek indices from the first portion of 
the alphabet a, 8, ... = 1, 2, ... The condition for phase equilibrium is 
simply the system of equations 



(18) 


Remember that we are suppressing the writing of the intensive thermodynamic 
variables which label the state. Let k be a fixed index and examine the 
behavior of dilute solutions in the vicinity of the point xj = s-jk, that 
is, solutions in which species k is dominant and all species j £ k are 
trace species. Then it is an experimental fact, known as Raoult's law, that 
the mole fraction of a trace species in one phase is proportional to its mole 
fraction in some other phase, assuming that we are sufficiently close to 
x-j = 6 -j k • Of course, the proportionality depends on the choice of trace 
species, dominant species, and the pair of phases being considered. 


a ^8 3 . , . 

*j - Vj J * k, as *, * s ik 

A natural generalization of this is the equation 
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(19) 



j ^ k, as x i ->■ « ik 


As the concentrations of the trace species increase, K 


iB 

jk 




and 


will 

J 


become concentration dependent. In the thermodynamics of multiphase systems 
one usually assumes that the phases are noninteracting in the sense that the 
Gibbs free energy of the whole system is the sum of the Gibbs free energies 

aB 

for each phase. This implies that one can factor into contributions 

from each phase. Thus we write the factorization 


nt = K*/K“ 
jk jk jk 


( 20 ) 


This factorization is not unique because we can put a common factor in the 
numerator and denominator on the right side and not affect the left side of 
equation (20). Let eju be any quantity with dimensions of energy per 
mole which depends on trie choice of dominant and trace species but not 
phase. Then by taking the logarithm of equation (19) and combining it with 
equations (18) and (20), we obtain 


e jk(<; + in K “ k ) - “j - e jk( q j in + ,n K jk) 


jkj 

j ^ k, as 


x . 


‘ik 


( 21 ) 


From equation (21) it is clear that we will obtain the condition (19) from 
the solution to the equations (18) if we equate both sides of equation (21) 
to a quanity fjk which depends only on the choice of dominant and trace 
species. 


- L* + e jk K j k * 


ct -| a 

e jk q j 


j ^ k, as 


x . 

i 


6 ik 


This form for the chemical potential can also be extended to the dominant spe- 
cies because as x^ + s.^, In x^ + 0 and the In x“ term does not contribute. 
Hence we tentatively write the chemical potential for all species in the form 

“j ■ “jk + 'jk x j as X i * «ik < 22 > 

where the phase index has now been suppressed and ujk and ejk may be 
composition dependent. 

Equation (22) is a nice simple form for the chemical potential as 
xi + 6-jk and is consistent with the generalization (19) of Raoult's law. 

But this is not sufficient. The chemical potential must also be the gradient 
of the Gibbs free energy, that is, there must exist a function g which 
generates equation (22) by differentiation as shown in equation (9). But 
this can only happen if the chemical potential satisfies the integrabi 1 ity 
conditions 9uj/3n|< = ank/anj (see, e.g., theorem III. 8, p. 1600, ref. la). 
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If we temporarily treat ujk and ejy as constants and apply the second 
member of equation (16) to (22), we oBtain 

n( 4 "j /in i> ' e jkV x j ‘ 'Jk 

Thus we find that 


n 




ik 


'jk 


( 23 ) 


and the form (22) does not satisfy the integrability conditions even with 
ujk and ejk as constants because the trace species i and j could be- 
have quite differently near the point x Si = 6 ^ and therefore we do not wish 
to assume that e ^ and ejk are the same. This forces us to look for the 
simplest modification of equation (22) which does satisfy the integrability 
conditions when yjk and ejk are independent of composition. Observe that 

9 N N 

n 7n7 ^ x «. e K.k = e ik ~ ^ x i, e a,k 

1 jL J. A» — J. 


and thus 


n 




x £ e £.k e ik 


'jk 


(24) 


Equations (23) and (24) show that if a linear combination of the e^k is 
subtracted from the right side of equation (22), a function results which 
satisfies the integrability conditions. Any quantity independent of mole 
fractions could be added without affecting integrability and, since adding 
e jk will give a simpler function for the Gibbs free energy, we write 


y j " p jk 


'jk jk 


In x. - £ 


a=l 


Vlk 


as x . 


ik 


(25) 


From this function the Gibbs free energy per mole in the vicinity of 
x-j = 6 j k , G(k), can be obtained by integrating the first member of equa- 
tion (16) or more simply by using the last equality in the first member of 
equation (14). Its validity can be checked by differentiation. 


G(k) - L (p jk * e jk 


x i )x j 


as x.j 


6 ik 


(26) 


Finally, to obtain a function which is suitable away from x-j = 6-jk, I shall 
merely superpose the G(k) weighted by the mole fractions. 


N 

G = 2 G(k) 
k=l 



N 

£ ( “jk + e jk ,n x j> x j x k 


(27) 
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The simple function (27), with yj^ and ejk independent of composi- 
tion but functions of the intensive variables, should represent the behavior 
of an N-ary system for the dilute solution regions in the vicinity of each of 
the N points Xi = 6-jk, k = 1, 2, ..., N. Consequently it is probable that 
it will also do quite well in other parts of the domain of the function. 
However, it is still possible to improve its flexibility and enhance its 
capability for dealing with concentrated solutions without diminishing its 
capabilities in the dilute solution regions of the domain. This can be ac- 
complished by the simple expedient of returning the parameters yj|< and ejk 
in equation (27) to composition dependent functions. We must then face the 
problem of making an explicit choice for the functions yjk and 
Clearly, one might concoct many potentially suitable functions on the basis 
of some selected criteria or personal whim. One can, in a sense, avoid mak- 
ing an explicit choice and permit the thermodynamic system to make the 
selection by working with an expansion of yjk and ejy as a linear com- 
bination of a suitable finite basis set of composition dependent functions. 
The expansion parameters can then be determined from data. That is the ap- 
proach I shall pursue here. I shall only demand that the character of the 
expansion permit the functions yjk and ejk to approach constant values 
in the dilute solution regions so J that conventional Raoult's law behavior can 
be accommodated there. Probably a number of basis sets exist that one might 
wish to consider. But there is one which comes to mind first because it 
offers both simplicity and elegance and at the same time is consistent with 
the six criteria enunciated at the beginning of this section: the set of N 
elementary symmetric functions of the N variables x-j. The elementary sym- 
metric functions <p^ £ )(x-j; N) are most conveniently defined by their genera- 
ting function ^(z; x-j, N) . 

N N m k 

#(z; X., N) s TT (1 + X,z) = E ?' (x . ; N)z (28) 

1 k=l K k=0 1 

The functions cp( i ) can be obtained from G by differentiations with 
respect to z. 

<p ( a) ( x . ; N) = &(z; x-, N) (29) 

1 a l dz* 1 z=0 


Clearly cp(0) = 1 is trivial and of no interest. For N = 3 the nontrivial 
functions are 


ce^^x.; 3) = x 1 + x 2 + x 3 


<p^(x^; 3) = x 1 x 2 + x 1 x 3 + x 2 x 3 V 


<p^ 3 ^ (x^ ; 3) = x 1 x 2 x 3 


( 30 ) 


The elementary symmetric functions of the mole fractions have some con- 
venient properties. Suppose we consider an M-ary subsystem of the N-ary 
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system. Then only M of the mole fractions are nonzero, and it is obvious 

( £ \ 

from the definition of the generating function that <p v '(x.; N) = 0 

for M < i <_ N. Furthermore for 1 <_ z <_ M the cp v y (x. ; N) are identical 
with the corresponding elementary symmetric functions formed from the M 
nonzero mole fractions. This can easily be verified from equation (28) or 
from equation (30) for the special case N = 3 by setting X 3 to zero, 
for example. A generating function for the partial derivatives of can 

be obtained by differentiating the generating function (28) with respect to the 
mole fractions. 


8 #(z; x^, N) 
_ 


k =0 


k 3 <p^(x.; N) 

3xT 


But 

9 In 9(z; x i , N) z 

3 X^ = (1 + x a z) 


(31) 


and therefore we obtain a very simple result. 

3 ^(z; x,, N) z^(z; x,, N) 

= = z TT (1 + x. z) = z^(z; x,, i 4 i t N 

3X * (1 + x £ z) kH k 1 


1 ) 

(32) 


If we apply the Liebnitz formula for the differentiation of a product uv in 
the form 


(uv) 


(n) 


£(> (k, * <n - k) 

k =0 VK/ 


to the function z^, we obtain 

(zSf)^ = + z0^ 

But then 


3«p^^(x i ; N) 
— 



i ± N 




z =0 




• (j - rrr x i> 1 * *• N - uj 


z=0 


and so we have shown that 
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(33) 


<pi j) (x i ; N) , 


3cp^ J ^(x,; N) / • n 

i 1 * i, N - 1) 


3X„ 


All higher derivatives can be obtained by further differentiation of equa- 
tion (33). The validity of equation (33) can be verified for N = 3 by 
direct differentiation of equation (30). 

Now that the properties of the basis set have been established, it is 
simple to write the final form for the function which will be used to repre- 
sent the Gibbs free energy in an N-ary system. We write 




► 


J 


(34) 


f A. ) f \ 

where vt\^' and are independent of mole fractions. Substitution of 

this expansion into equation (27) gives the appropriate function for the 
Gibbs free energy in expanded form. 


G = 







x jX k 


= £ <P U) (x.; N)G U) ( Xi ) (35) 


This representation for the Gibbs free energy defines a new class of non- 
ideal solutions whose properties will be studied for the balance of the 
report. The structure of equation (35) shows that only the symmetric part 


of 



can contribute to G, and so without loss in generality we take 




(36) 


This property does not apply to e 


(i) 


because the In x. 


. • * ucLQusc unc iii a • term introduces 

<J J 

asymmetry. If the full expansion (35) is truncated after the first L 
elementary symmetric functions, then the expansion will be called an expan- 
sion of degree- L. The number of independent parameters for each «. is 
N(N + 1)12 + N^ = N(3N + 1 ) /2 , and hence for an expansion of degree L there 
are LN(3N + 1 ) / 2 independent parameters. 

The properties of the elementary symmetric functions allow a simple in- 
terpretation of the terms in the Gibbs free energy expansion (35), which 
correspond to different values of a. The functions <?(*■) can be regarded 
as weighting functions for different regions of the domain of 
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N (1) ^ 

definition 0 < x. < 1, E x. = 1. Since 9 V = E x. 

(1) k=l k 1=1 

9 V ' is a constant function over the domain and thus the 


= 1, the function 
a = 1 term 


samples all points of an N-ary system with equal weight. The function 

9 ^ = E x -x . vanishes at the N points x - ( k ) = s.. , k = 1,2, N, 
i< j J i ik 

and so the a = 2 term samples all points with the exception of those N 
points corresponding to the unary subsystems. The function 9 ( 3 ) vanishes 
for all unary and binary subsystems but samples all other points. Similarly 
vanishes for all unary, binary, and ternary subsystems. In general, 
r 4 )(xj N) is zero for all M-ary subsystems for which M < a. This 
situation can be visualized by the description that, as a increases, the 
terms progressively sample more and more of the "interior" and less and less 
of the "boundary" of an N-ary system. This description is especially appro- 
priate when applied to the graphic portrayal of ternary systems by equilat- 
eral triangular diagrams where the vertices are the unary subsystems, the 
sides are the binary subsystems, and the interior is the portion of the ter- 
nary system for which (x-j ^ 0, i = 1,2,3). 

There still remains the task of calculating the chemical potentials for 
the Gibbs free energy function (35) by using the last equality in the first 
member of equation (16). Suppose we define 


.(*) = 


(1 


0 . )G 

i ' 


(1) 


(37) 


were the differentiation operator 0 ^ is defined in equation (16). Then 
from the last equality in equation (35) we find 


V 1 J (*•)/■ .1 \ (a) , r (a) (a), v 1 (*>)/ 

v i " & 1* < x j ; N >“i + G *1 (x j ; N > - Vk < x j ; 


N) 


* f £, ) 

But since 9 Ms homogeneous of degree a, we know that 
N / \ 

T f = S -9 ^ and thus 


k=l 


u i = E j^J - aG^j 9^(x^; N) + G^ £ ^9 ^(xj; N) 

(1) A 

The i=l term can be simplified because 9 V = E x . = 1 and 

j =1 J 


39^/8X.j = 1. 

Thus we can write 
N 


- JV + 


y i = 


E 

1=2 


(u ^ - aG^^ 9^^(Xj-; N) + G^^^^Xjj N) 


(38) 


(39) 
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By using the definition of y\ ' (eq. (37)) and the symmetry property 

( i) 1 

(eq. (36)) of the parameters yi\. ' it is a simple calculation to establish 


an explicit formula for 

4‘ > - s[ 2 " ( ik’ * s 1k’ (1 * ln x i> + 4?’ ln x k]> 


'jk 


■ 5 k?i[^ + (1 + ln Xj) ] XjX| 


■ ( £ 1 [ 2 4k ) + 4-k 1 (1 + in x i> + 4i’ in x J x t 


- G 


(l) - E E 4 k l> Vk 

j=l k=l JK J K 


(40) 


It is sometimes convenient to have an expression for 
grouped by parameter. 


.(*> < 


in which terms are 


N N f 

= E E i x i (26 ik " x k )y 

1 j=l k=l ^ J 


U) 

jk 


(x j 6 ik + ViJ - x j x k } ln x j + x k< 6 ij 


X j> 


.(*) 

'jk 


(41) 


At first glance the result given by equations (39) and (41) is a bit discon- 
certing for it seems to imply that we can extract all parameters in the Gibbs 
free energy by fitting one of the y-j alone. This seems to contradict the 
fact that differentiation of a function results in a loss of information 
about the function because in the differentiation we discard the constant of 
integration. That, in fact, is_also the case here. To see this most simply, 
we need only define functions y-jj(x k ) by 

"i.i (x k> - ?i.i< x k> - <"i * 


and the n k are arbitrary and independent of composition. If this is sub- 
stituted into the expression for the Gibbs free energy (eq. (27)), we find 

N N N N _ N 

£ £ <u J k + ‘j" ln x J )x j Xk ■ £ k£ ( “ 3k + 1n Xj)XjXk + £ 'ft 


If is the chemical potential calculated from the yj k and 


is the 


chemical potential calculated from the Uj k , then the two are obviously re- 
lated by the expression 
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u 


i 



From this it follows that there are N - m arbitrary parameters in the Gibbs 
free energy function (eq. (27) or (35)) if we have data on only m chemical 
potentials. 

It is occasionally of some interest to have an expression for the sec- 
ond derivatives of the Gibbs free energy. An appropriate formula for the 
calculation of the second derivatives from G can be obtained easily from 
equation (16). 


n 



Vi - V 1 + V G 


(42) 


The left side is manifestly symmetric, but the symmetry of the right side is 
not obvious. To demonstrate the symmetry in a simple way, it is convenient 
to use the commutation properties of the operators ^ k . One can evaluate 
the commutators of the operators from the definition of the operators by a 
straightforward computation. 


rW * 




- = 
i k 


8X. 



C*. 


(43) 


Symmetry can now be simply demonstrated by combining equations (42) and (43). 


n 




c/yi + 


+ ^ k )]G 


= (^ k - C s + O k , ^])G = 0 

The actual calculation of the second derivative by means of the first equal- 
ity in equation (42) can be somewhat simplified by taking advantage of some 
properties of the operators /? k . When /? k is applied to a product of 
functions, it has a simple behavior. 

^ k (uv) = u *? k v + v^ k u (44) 


When it acts on the functions G^, q>^, and <p^, we obtain 


€ 

K IN 

0 U) (a) „ (a) 

e 'k <p = <p k “ ^ 


«,(*) - u - d,,! 1 * 

k Y i 3x k ax^ v i 


> 


(45) 
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The first of these three results is an immediate consequence of equation (37), 
and the remaining two follow from the definition of the operator and the fact 

that 9 ^ is homogeneous of degree The effect of C ^ on can 

be calculated by expressing ^ ^ in terms of derivatives of G^. 


Vi*’ - ‘V 1 * ^-) G 


( a) 3 G 


2 C (M 


3 3X 


riM 


A^> 

Sx k SXj 


Ail 

3 X • 3X. 
i J/ 


N N 

* Z Z v 


a 2 G< l > 


J m 3X. 3 X 

j=l m=l J m 

The evaluation of the right side of this expression gives the final formula 
which is required for the evaluation of the second derivatives of g. 




k y i 


(0 


o V *• 
= . 


kf * 44 d + in x k> + 44 i 1 * in x i> * x 4 6 ik E 44 x j 

vJ ^ 


- £ - 44) * 45’ f 1 + in x k) + ‘'iZ (i * >» x i> 


(4k 1 + 44) ♦ 1- x :)] * ^ | i [44 


l y mj + e ij } t 1 + ln x m) 


)] XjX m 


(46) 


The second derivatives are obtained by combining equations (42), (44), and 
(45) with the expression for y-j given in equation (39). 


a 2 g 


an^ an 


r W | 2 {* (l) V l) - H ,, + 4‘ , )(*’“ , -4‘ ) ) 


44 Z-_ ( , + _ < t . d^ 1 ) + m ♦ D» 


(*) 


:(«-) 


(47) 


The function I have constructed to describe the composition dependence 
of the Gibbs free energy for a class of nonideal solutions contains an impor- 
tant special case. Suppose that we make a special choice for the functions 
yj k and ej k , which appear in equation (27). 


p jk = u jk E (u j + v k ]l2 


( 48 ) 


e jk " e jk 
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Then 


5 5 | k.l ( "JK + 'J k ^ X J >X j X k ‘ j| ( "j + * x j> 


,)Xj 


(49) 


and if e and the y. are taken to be independent of composition, then 5 

J ^ — 

becomes the ideal Gibbs free energy if e = RT. The function G is useful 
for defining a relative Gibbs free energy function which has the same form 
as G itself. 


N N 

aG s G - G = £ 2C (Ay.j k + A Ej . k In Xj)x j .x (< 


N 


Au jk = 


Ae jk " 'jk - 6 ' £ “jk )<p< ‘ l|(x i ; "> 


* “jk - “jk * “jk - ( “j + “k > /2 * E N > 


4=1 


> 


(50) 


Naturally aG becomes identical to the excess Gibbs free energy aG 
e = RT and 


if 


u . — y • * 

j j 


where y. are the chemical potentials of the pure con- 
J 

stituents, that is, the Gibbs free energies per mole. This means that the 
expression for the derivatives of the Gibbs free energy are immediately 
applicable to the excess Gibbs free energy merely by replacing y - k and e- k 


with 


Ay jk 


and 


Ae j k with the understanding that 


r 


Ay 


(0 

jk 


= y 


(a) 

jk 


- (^j 


y k 


)/2 


(a) 

= y 

v Jk 


4 = 1 


l + 1 


»_ ( ^ ) _ ( ^ ) DT 

A jk - e jk _ RT 


The special choice Ay:; 


= 0 shows that 


jj k ^ =0 for 4*1 and Aej^ 

Hildebrand's regular solutions are special cases of the class of solutions 
considered here. 


GIBBS FREE ENERGY IN A SUBSPACE 


I have already briefly mentioned M-ary subsystems of an N-ary system 
(M ^ N) in the preceding section. This is nothing more than a special case 
of another more general notion, that of a subspace, which also is relevant to 
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thermodynamics. The simplest example is an aqueous sodium chloride solution, 
which is a ternary system with constituents H 2 O, Na + , and Cl. The solution 
is formed by dissolving NaCl in water and therefore the moles of Na + and Cl- 
are always equal. This implies that there are not three, but only two inde- 
pendent composition variables despite the fact that we have a ternary system. 
The independent variables are clearly the moles of water and the moles of 
sodium chloride. The form of the Gibbs free energy written in terms of the 
independent variables will differ from the form written in terms of variables 
which are not all independent. In this section I shall examine the form of 
the Gibbs free energy in a subspace of an N-ary system. 

To carry out this task expeditiously, I must introduce some helpful 
notation. Latin indices will continue to use the range 1, 2, ..., N, while 
Greek indices from the last part of the alphabet (u , v, a, t, ...) will use 
the range 1, 2, ..., M £ N. For thermodynamic purposes an M-dimensional sub- 
space can be defined by the equations 


n i - n i< n 0 > 

n .( An ) = xn . ( n 
v a ' i ' c 



(51) 


where, of course, n a represents the independent composition variables of 
the subspace in moles per unit mass and the second member of equation (51) 
is the requirement that the functions n-j(n a ) be homogeneous of degree 1. 
The functions n-j(n 0 ) always turn out to be linear functions in actual 
situations, but this specialization does not lead to any mathematical sim- 
plification, and so I shall use the more general case. Suppose we define 
"stoichiometric coefficients" v-j a by 


an . 


'-I, ■ ”i 0 (n ,) s Irf — 0 


(52) 


where v ia would be a constant in the usual case of linear functions. We 
can now apply Euler's theorem to equation (51), as we did before in the case 
of the Gibbs free energy, to obtain 


n. = £ v. n 
1 1 a a 

o=i 

v- (xn ) = v. (n ) 

To' T ; 1 0 ' T 




(53) 


The total moles n and mole fractions x of the independent variables are 
defined analogously to equation (13). 


n 



n 


a 


> 


X = 

T 


n /n 
x 




(54) 
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If we define 


V 

a 



> 0 


(55) 


then we have simple relationships between the two kinds of total moles and 
mole fractions. 


x 


i 


- -1 
nn 


M 

y v, x 

M loo 
o=l 




v x 

T X 





(56) 


We are now in a position to define the Gibbs free energy functions and 
the chemical potentials in the subspace. 


9 = 9(n o ) = g[n i (n a )] ^ 

r A -1 A / \ 

G E n g(n 0 ) 



(57) 


The function g is homogeneous of degree 1 in the n CT because g is homo- 
geneous of degree 1 in the n-j and the nj are homogeneous o£ degree 1 in 
the n a . It is simple to establish the relationship between G and G and 
between y CT and p j . 


G = 


n 1 g[n i (n c 


) ](n/n) = Gn/n 


u 


a 




M _-v. 

1 1 a 


(58) 


A complete parallel exists between the N-ary equations and those of the sub- 
space. The counterparts of equations (14), (16), (43), and (44) are 
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M 

G(x p ) = E Va 

a=l 


% = ^ - E X T tt; + Trj 2(x p } = (1 + ^(* p ) 


3u 

a 

an 

V 



^ P (x ) 
v o' p' 


^ (1 + 0 )G(x ) 

V v' cr ' p' 


( 59 ) 




^1- 


^x- 


^ CT (uv) = u ^? a v + V /? a u J 

While the expressions given in equation (59) for y a and 3y a /3n v are 
correct, they are not necessarily the simplest ways to calculate these quan- 
tities. It might be simpler to calculate y a from the second equation in 
(58) if one already has expressions for y-j. A similar formula can easily be 
found to relate the derivatives 3y a /3n T to Sy^/snj merely by differentia- 
ting the second equation in (58). 


3y 

CT 


an 


T 


N N 

v v,.,. 

Z—f Z_-/ JT 1CT 

j=l i=l 


3y . v 3v . 

+ y . — 

3n . 1 3n 

J i = l 


(60) 


Only the first part of equation (60) will contribute in actual problems be- 
cause, as I pointed out before, 3v-j CT /sn T = 0 in the usual case. The Gibbs- 

M 

Duhem equation for the subspace, E n 3 p /3n =0, follows from equations 

a=l 

(60) and (53) and the Gibbs-Duhem equation (12). 

From equations (35) and (58) it may be seen that G can be written as 

G = E^ «p (£ ' ) [x.(x o ) ; N]G U) (61) 

where 
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8<1> 5 Sn_1 5 li i*" * J *j_ 

(a) _ 'f v' (a) (a) 

w' 7 = 2^ 2- v. v. u' y = ui ' 


X 0 X T 


>1 kti J« k ^« " 


,(*> = y v ,(0 

' jT ' k=i kT jk 


(62) 


Similarly the combination of equations (39) and (58) gives an expression for 
^0* 

(1) . f( (a) .<* -l~(i.)\ ( a ) (- / . -li(i) ( a)1 

v 0 = u CT + 2J 1^% - v a Ann G v y y [x i (x £j ); N] + nn G v 'q>‘ 'j 


(*■) . y U ) 

'« = fi W i V - 


> (63) 


,(*) 


N 


fi 


U) 


V ia 




To obtain an explicit expression for j/ , we can use equation (58) with 

either the first equality in (40) or (41). I shall only give the one which 
comes from equation (41) since it gives the most compact form directly. 


v 


(*) = 


r» n i 

E ESrf 1 * 

P =1 T=1 



A — 

nn 


*v X ) 
0 t/ 


y 


U) 

PT 



In 


x i + 


- -1 , 

nn x (v. 

t' lo 


X 


,X1 



G 


U) 

i t 


M M 

= E E 

p=l T=1 


* -1 
nn 


(» 


6 + 
P 0T 


X 6 

T 0p 



v X X 
0 p 



N 

E 

i=l 


m r- 

E P 

T=1 L 


s x . 

GT 1 


In 


x i 


nn 


-1 


x ( v . - 

T ' 10 


x i v o ) ( 1 


+ In 


il (*) 

X i } J £ iT 


(64) 


The second equality is obtained from the first by using the symmetry of the 
(eq. (62)). The number of independent parameters in G is easily cal- 

pt ( i) ' 

culated since for each a we have M(M + 1 ) / 2 independent u v ' and at most 
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MN independent e^. ^ and hence for L terms the total is LM(M + 1 + 2N)/2. 
1T ( 

The number of independent e- ' parameters might be less than MN for each 


i for special choices of the stoichiometric coefficients, 
a given i and j we have v. = v. for all 0 . Then 
therefore for all t 1ct JtT 


Suppose that for 
x • = x • and 

* J 


x 

T 


M 

£ "i/ 

o=i 


1 n x . = x 

1 T 



v. x 
jo a 


In x. 

J 


( 

and therefore the functions multiplying e^ T 
which precludes the separate determination of 
case, for all t and only the combination 


, (1) 

and 

^ and 

0 


in 


.( M 


IT JT 

( 0 > ♦ ■$:>) « 


can 


are equal. 
In such a 
be deter- 


mined by fitting data, and one may impose arbitrary subsidiary conditions to 

determine separately the individual coefficients. 

One can introduce a relative Gibbs free energy in the subspace in much 
the same way as was done for an N-ary system by equation (50). 


* a _i 

aG = nn 


M M 


J2 [am + Ae. v. In x. X X 

0=1 ti \ OT jti Ja 3 0 T 


N 




- MO U), 


5 0 + n x )/2“ £ Au-\-^[x i( x p ); N] 

5*— J. 


(o5) 




N] 


The relative Gibbs free energy aG becomes an excess Gibbs free energy if 

e = RT and y = y* where y* are the Gibbs free energies per mole for 
0 o’ 0 3 

the independent species. I should point out that Ay pi . and Ae^ are n °t 

the same as the transformed versions of Ay ., and Ae ., , which are given in 
equation (66). 3 3 


N N 

E £ S v ja v kT^jk = 

N 

Ae . = /. \i, Ae = e. - ev 

JT kT Jk JT T 


-(yv +yv)/2 
v t H t a 


N 

V = S v ,- VI. 

jTi J 


( 66 ) 
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The differences between am , Ae. and Ay , Ae. can be seen by 

OT* JT O T» JT 17 

comparing the last two members of equation (65) with the first two members of 
equation (66). 

The choice of independent composition variables is not unique and often- 
times one choice mjght be more convenient than another. Suppose we wished to 
use the variables n 0 in place of n 0 . Then for thermodynamic purposes the 
transformation between the two sets of variables must satisfy 


n ( An ) = An ( n ) 
cr' p cr ' q 



where J is the Jacobian of the transformation. It is an immediate conse- 
quence of the chain rule for differentiation that 


an i 


an 


a 



v. T 
ip pp 



M 

= E 

P =i 


u T 

P PO 


and obviously 


;r<‘> 

QT 


M 

= E 

p=i 


M 

E 

A=1 


y( A )j T A 
pA per At 


with similar transf ormation rules for other quantities. I shall not pursue 
this farther except to mention that the T pcr will generally be constants. 

LIMITING BEHAVIOR 

Thermodynamicists are often preoccupied with the limiting behavior of 
thermodynamic properties in dilute solutions where the concentration of one 
species dominates that of all other species. This naturally raises the 
question of the limiting behavior of the Gibbs free energy function and its 
derivatives for the class of solutions discussed in the preceding section. 
For this reason I shall examine the dilute solution behavior of these solu- 
tions and shall obtain an interpretation of some of the parameters as a 
byproduct of the analysis. The limiting process will be carried out for 
constant values of all intensive parameters, such as temperature and pres- 
sure for fluids, which characterize the thermodynamic state. It should be 
kept in mind that although this limiting process always makes mathematical 


k 
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sense, it might not always be physically meaningful. For example, suppose 
that we ask for the limiting behavior of an aqueous sodium chloride solution 
as the concentration of NaCI goes to zero at some given temperature and a 
pressure less than the vapor pressure of pure liquid water at that tempera- 
ture. But pure liquid water is inaccessible from the aqueous solution phase 
under the specified conditions, and the limiting process is not physically 
meaningful. Similarly the limiting process is not physically meaningful in 
the pure sodium chloride limit at ordinary temperatures and pressures 
because NaCI has a limited solubility in water. The formulas which follow 
all presuppose that the pure state x-j = x-j(p) = 6-j p , where p is a fixed 
index, is physically accessible from the phase under consideration and, 
strictly speaking, they are physically meaningless if it is not. 

To simplify the calculation of the limiting values of the Gibbs free 
energy and its derivatives, it will be helpful to know some auxiliary 
limiting values. 

■N 

1 im x a In x = 0 for all a 
x->-l 


lim x In x = 0 a>0 

x-*-0 


1 im 

x r 6 i P 


x j 6 jp 


i N > - 5 u 

* (1 


ax. 


s kp> s 2* 


v. 


> 


j 


(67) 


The second of these five limiting values comes from a simple application of 
L 'Hospital's rule and shows that In x has a weak singularity indeed in the 
limit x -*• 0. The fourth and fifth values make explicit use of specific 
properties of <p(l) and <*>(2), namely. 



N 


E x. = 1 
i=l 1 


9 


( 2 ) 


E 

i<j 


X i X j 



1 


> 


( 68 ) 
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and that <?(*•) is of degree ft in the mole fractions. We now use the 
limiting behavior given in equation (67) in evaluating the limiting behavior 
of G as given in equation (35). 


t * = lim G = lim 


:(1) _ ..(1) 


= v 


V 5 ip 


*i* 5 ip 


pp 


(69) 


Hence we have the important conclusion that when a pure species is accessible 
from a phase, the corresponding diagonal parameter must be chosen as the 

Gibbs free energy of the pure species. From the expression (41) for p^ 
it is easy to see that 


lim 

x -*-s 
m mp 


,0 0 


N 

E 

k=l 


(ft) , 
e\, 'x. In X . 
ik k i 


= lim 
x 

m mp 


.(*) 


N N 

E E 

j=l k=l 


(a) 

e jk Vij ln x j 


N N 

= E E , 

j=l k=l L 


( ft) 

6 jp^ 26 i k " 6 kp^jk 


+ V 5 U - 


U) + (*) . <«) + .(»>_ ,(») 

^1P M P1 M PP IP PP 


(70) 


N 


where, of course, in ]£] e^x. 1 n x. the term k = i vanishes in the limit 

k=l lk K 1 


as does the complete sum for i = p. Suppose we define 

N / . x N 


w,(p) = lim 

X +6 

m mp l 


w, - E ■P U) E ^>x k In x, 

1 ft=l k=l lk k \ 


then it follows from equation (39) that 

N 


n,(p) = lim 

1 x +6 

m mp *- 


(1) V' n . r (2) (2) 
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and from this we see that for i = p 
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The chemical potentials in a dilute solution have a simple form as can be 
seen from equations (71) and (34). 


Ui = u“(p) + 
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£ 

k=l 


G ik x k 


In x. 


as 
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(74) 


If we compare this to the defining equation for the rational osmotic coef- 
ficient (Kirkwood and Oppenheim, ref. 23, p. 164, eq. (11-31), and Pitzer 
and Brewer, ref. 24, p. 321, eq. (22-22)), we see that the similarity of 
form permits us to interpret the functions e-jj as generalizations of the 
rational osmotic coefficients. 

The limiting behavior in a subspace is a bit more difficult to discuss 
because we must now take into account the stoichiometric coefficients v-j p 
and because we are now concerned with the behavior as x p + x p ((o) = 6 , 

where <o is a fixed index, and not as x-t + 6-j p . Instead of the relatively 
simple limiting values shown in equation (67) we have 
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where in equation (75) and for the remainder of this section the symbols v- 
are to be interpreted as lim v. . The second and third parts of equation 

x +6 u 

p P (0 

(75) are a direct consequence of equation (56). The fourth part can be ob- 
tained simply by looking at the limiting behavior of the generating function 
^(eq. (28)) and recalling that the q>( a ) are homogeneous of degree a, while 
the fifth part comes from equation (33). The last part of equation (75) is a 
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consequence of Euler's theorem applied to From the definition of G^) 

given in the first member of equation (62) we see that 


lim G^ 
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and so from equation (61) 
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and as a general rule y* ^ y(^) although it could still be true in some 

particular case. One very important special case for which it is true is 
the one corresponding to an aqueous electrolyte solution in which water is 
assumed to be nonionized and the index w labels water. Under these condi- 
tions the stoichiometric coefficients for water satisfy v. =0 or v. = 1 

JW JW 

( 1 ) 

and v = 1. For these stoichiometric coefficients only y v survives on 
w ww 

the right side of equation (76). In general, from the nonnegativity condi- 
tions (52) and (55), it is sufficient to require only that v u = 1. 
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★ 

0) 



if v = 1 
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The condition in equation (77) would not be satisfied if we took the ioniza- 
tion of water into account. 

Next I shall establish the limiting form of the chemical potential y 0 
in the vicinity of the point x p (w). The chemical potential y j only had 
the logarithmic singularity In x-j ; however, y CT can have a logarithmic sin- 
gularity for each i for which v-j a = 0. Therefore, in contrast to equa- 
tion (71), I shall define 
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From the first expression in equation (64) for y^ we have 
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where I have used notation analogous to that in equation (62) 
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If we use the expression for given in equation (63) and take advantage 

of the fact that <p^ = 1, then 
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From this equation and (79) it immediately follows that 
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and this result is considerably more complex than the corresponding formula 
for y ”(p) given in equation (72). From equation (78) it follows that 
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0 
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52 52 v . e . x In x . 

f-s 10 IT T 1 
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represents the limiting behavior of p a . 

The expression (80) for p°°((o) simplifies considerably for o = <d. For 
this case the second summation 8ver 4 vanishes identically by virtue of the 
last member of equation (75) and only the first summation over 4 survives. 
From this we see immediately that 
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Thus, in general, y”(<o) ^ p* and from equations (82) and (76) we can obtain 
an explicit expression for the difference 
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This formula can also be obtained directly from equation (81) with a = u 
and u* = lim p . It is clear from equation (83) that p°°(d>) and p* 

u „ (0 d) ' d) 

X +5 
P Pd) 

will be equal only under special circumstances. One such special case is the 
condition v = 1, which implies that v. /v =0 or 1 for all j values, 

0) Jd) d) , 

but it could also be true under other circumstances since the e\ ' may be 

jd) J 

positive or negative. The assumption = 1 also engenders a considerable 

simplification in the expression for p”(o)), o ^ to. Under these 
circumstances 
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These results, when substituted into the formula (80) for produce 

the expression c 



2 U (D + ^ 

aw aw 




a 4 = <d. 
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1 


(84) 

where the sum over k contains only one nonzero term because v =1. Of 

course, if the stoichiometric coefficients v. and v. , a ^ u>“ are 
’ ku Ka 

orthogonal, as is generally the case, then the sum over k vanishes. 

Thermodynamic properties of solutions are always reported relative to 

reference values which, of course, may be chosen arbitrarily. Thus one would 


express the chemical potentials for an N-ary system as (y.. - y.) or as 
(u - u ) in a subspace, where y. and £ are independent of composition. 

a a la' 

For the class of solutions described by the equations of this report one 
might reasonably select = y“(p) and y^ = y”(w). Perhaps the most com- 
monly used reference values, especially for electrolyte solutions, are the 
so-called infinite dilution standard states. These standard states are 
described, for example, by Kirkwood and Oppenheim (ref. 23). The use of 
such standard states is predicated on the existence of certain limits in- 
volving the chemical potentials themselves. The particular values of these 
limits are of no practical consequence, but their existence certainly is. 
The limits in question are defined by 


u°(p) = lim [y. - RT In x.] 

x -*-6 
m mp 

u^(w) = lim 
x ■+& 

P pin 



(85) 


( 86 ) 


where m T are molalities relative to species w. The first of these is 
Kirkwood and Oppenheim's equation (11-8), while the second is a combination 
of their equations (12-3) and (12-7). The infinite dilution reference values 

are y?(p) for an N-ary system and y*, y°(o>) for a 4= w in a subspace. Are 

the infinite dilution values suitable as reference values? Only if the lim- 
its (85) and (86) exist. These limits, like all limits, are experimentally 
unobservable and their existence for real solutions must be taken as an arti- 
cle of faith when one makes a decision to use them. Because the limits are 
not experimentally measurable, their values for real solutions can only be 
obtained by an extrapolation of data. But this makes the resulting values an 
artifact of the extrapolation technique and hence they can be of no intrinsic 
significance. Nonetheless, for the class of solutions defined by the equa- 
tions of this report and for those real solutions which can be adequately 
represented by these equations, it is possible to evaluate these limits 
exactly with the aid of the limiting behavior expressed by equations (74) 
and (81). 
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From equations (85), (74), and (67) we find that 
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This then permits the enunciation of a simple theorem: 


The limit p?(p) exists and y?(p) = u”(p) if and only if 
i = p or^eT^/RT = 1. 


The calculation of u°(<o) is a bit more involved. If we use m = x m /x , 

CT T T (0 (0 9 

then from equations (86), (81), and (75) it can be shown that 
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From this formula there results a second theorem: 


The limit u°(u) exists and y 0 (w) = y^to) - RT In 


v e /v RT 

(m ) 0 (v ) 

v or 'or 


if and only if v. =0 or (v. =0 and e. /v RT = 1) for all i 

lo ' 1(1) U) U) y 


For c = w we see that y°(u) can exist only if v. =0 for all i, which 

0) 1 0) 

clearly is not possible for it violates the condition (55) that v >0. 

Thus y°((o) never exists and y°(w), a 4 u, only exists under special circum- 

W O’ 

stances. The nonexistence of y°(a>) is a consequence of the composition 

10 

dependence used in the second term of the defining equation (86) and not be- 
cause of the behavior of y . The limits y°-(p) and y°(u>) generally do not 

w fsil 1 a 

exist simply because the parameters . are permitted to differ from RT 

for the class of solutions under discussion. 

Thermodynamici sts are not always content with giving thermodynamic 
properties relative to reference values. They often insist on expressing 
properties relative to a reference function with the expectation that 
properties relative to a suitably chosen function will exhibit simpler 
behavior than the property itself. This practice leads to the definition of 
activity coefficients both for N-ary solutions and for subspaces of these 
solutions. A common definition for the activity coefficient yj in an 
N-ary solution, relative to y-j, is 


m. s y . + RT In ( y ^ x i ) 


( 88 ) 


while for solutes in a subspace of an N-ary solution the activity coefficient 
y a , a 4 a), relative to y a is conventionally defined as 


y„ = 
a a 
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The following theorem is a trivial consequence of the definitions of the 
activity coefficients and the values y°(p), y°(u): 


The limit lim In y.- exists and equals y!?(p) - y i /RT if 

V -*.£ L * * J 


and 


m mp 


only if y?(p) exists. The limit lim In y , a 4= w, exists and 

1 x -*-6 0 

p pm 

equals (y>) - y g )/ v^RT if and only if p°(u) exists. 


Clearly, since y?(p) 


and y°(«j) do not generally exist, y^ and y a will 
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generally not have finite limiting values for members of the class of solu- 
tions under discussion. Nonetheless it is possible to define activity coef- 
ficients y. and Y a which will have finite limiting values for all members 
of the class. An obvious choice would be the definitions 


5 “>> + £ 'ik x k ln x i * RT ln n 


i N M 

y“(o)) + nn 1 Y\ v . 7] e..x. In x. + v RT ln T 

a la it t 1 a ’a 


(90) 


and it is clear from equations (74) and (81) that y. and y both approach 
unity for the infinitely dilute solution. 1 ° 

The result that y-j and y a may not have finite limiting values is of 
considerable importance because they are sometimes used to calculate practi- 
cal osmotic coefficients for the solvent in dilute binary solutions by quad- 
rature. The integration extends from zero solute molality to some finite 
value and requires the extrapolation of the solute activity coefficient to 
zero molality. Consequently the value of the osmotic coefficient calculated 
in this manner is dependent not only on the activity coefficient data, but 
also on the chosen method of extrapolation. The extrapolation methods chosen 
always assume that the solute activity coefficient (89) is finite, actually 
unity, at zero molality. The resulting formula for the osmotic coefficient 
(ref. 25) (also see, e.g., Pitzer and Brewer, ref. 24, eq. (23-22), p. 339) 
implies that its value is also unity at zero molality. But since the assump- 
tion on which it is based is not generally valid for the class of solutions 
under discussion, it is necessary to look directly at the limiting behavior 
of the practical osmotic coefficient for these solutions. 

Practical osmotic coefficients are defined in terms of the solvent chem- 
ical potential. For an N-ary solution we have the definition 

*p E - x p (u p " y 5 )/RT x i (91) 


while for a subspace it is 
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(see Pitzer and Brewer (ref. 24), eq. 20-28), p 
p. 570). If we use equations (91), (73), (74), 


263, and eq. (34-33), 
and (67), we see that 
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If we use the expansion ln(l + a) = a, it follows that 
^ 1 ^ / RT 
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To calculate the limiting value of we first observe that 
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But if equations (83) and (87) are used for p* - p”(<o) and equation (81) 
and the first member of (56) are used for p^ - p”(w), then we may write 


[p* - p (io)l - [p - p (J = -e v ^ In 
|_ O) to J |_ (0 (0 'J (0(0 to 


i A 

/ -1 > (OT 

v ( nn ) 
to' ' 


N 

- E 


°° -1 

e . v v . In 

-1 
V . 

/ M \ A jw 

( E V. X ) 

Jw to Jto 

Jw 

\ 0=1 JCT 7 


where 


M 


A i v nn ^ T] e x / e~ 

Cl) X 0) WT T (00) 

T=1 


as x +6 

P pw 


i 11 

A.. 2 v nn 1 E 


JlO (0 


T=1 


e _• X / e . 
JT T J(0 


But it is easy to see that 
lim A = 1 = lim A . 

10T 1(0 

X *6 X ■*'(5 J 

p pto p pu 

and consequently 


N 


— 1 ' oo 

lim <t> = -(v RT) E e ” v - lim 

x ->-($ 

P pw 


j=l X -6 , , 

p pto 


?E v x V 1 In (v' 1 E v.. A 
T 7 V a=l Ja 7 


If, for a given value of j, e. v. =0, then we qet no contribution to the 

Jto Jto ’ 

limit from this term. If, on the other hand, v. e . ^ 0, then this implies 

jto jo) 

that the term will contribute and also that v. ^ 0. For such a value of 
j we see, by using x = 1 - >, x , that 
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where in the last step I used 1 n ( 1 + «) s a. From this it is clear that the 
contribution of a particular term is dependent on the direction in which one 
approaches the infinitely dilute solution and consequently so too is the lim- 
iting value of <f w . Suppose we approach the infinitely dilute solution along 
a line of fixed solute proportions, that is, along a curve whose parametric 
representation in the vicinity of the infinitely dilute state is given, to 
first order in t, by 


x x = a xt, a^ ^ 0, S ^ ® 


where a^ are constants and t is a marking parameter for the curve. 
Clearly, we then have 
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If we designate the limiting value of <t> along this direction by <t>”(a, ), 
then u u x 
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(94) 


where the summation over j is to be read as the sum over all j values 
such that 4 0. For the special case of a path along which only one sol- 
ute is present we have aw^ = 0 and a a 4 0 if o is the label for that 
solute. Under these conditions the result (94) specializes to 


4> (a , = 0) = <(> (a) = (v v RT) 

to \4o a) o ti 
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If v =1, then only one term in the sum over j in equations (94) and (95) 

“ a, m 111 

survives and also e. = e\ 1 and e = e v From equation (95) we can 

Jo) JO) 0)0) 0)0) 

conclude that <f>°°(a) will depend on the solute only through v if and only 
if the sum over j vanishes. This will occur, for example, if = 0 
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or e“ u = 0 for all i. Similarly 
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CONCLUDING REMARKS 


In this report I have constructed a function to represent the composi- 
tion dependence of the Gibbs free energy for a class of nonideal solutions 
and studied some of its properties. The function has several useful charac- 
teristics: (1) Its parameters, which depend on the intensive thermodynamic 

variables which specify the state, appear linearly. (2) It contains a loga- 
rithmic singularity in the region of dilute solutions and contains ideal 
solutions and regular solutions as special cases. (3) It is applicable to 
N-ary systems and reduces to M-ary systems (M <_ N) in a form invariant man- 
ner. Because the parameters in the function occur linearly, all discussions 
of the Gibbs free energy can be translated immediately to other thermodynamic 
functions. For example, in the case of a fluid, the intensive thermodynamic 
variables are the temperature T and pressure P, and the Gibbs free energy 
(eq. (50)) can be converted into equations for enthalpy and volume by 
differentiation. 
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The utility of the function I have constructed must be judged not only by 
its convenient analytical properties, but also by its ability to reproduce 
the data for difficult real systems. This aspect of the problem is examined 
in a companion report (NASA TP-1930), where we determine the parameters from 
experimental data for several systems by the method of least squares. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, December 15, 1981 
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